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ABSTRACT 

peewo-level numéeéuvediepwcavetron model incorporating 
terrain and non-adiabatic warming effects is tested. An 
attempt is made to determine the empirical coefficient (s) 
EO~, Cach of three terms compesing a prognostic Gquation 
designed to predice eae 1000 toms0C-meseaitekness Ehrough 
Ghe use Of a digital ecompucer. 

The model is applied to three cases in April for a 
large part of the Northern Hemisphere. The predicted 
positions of most of the pressure systems were reasonably 
accurate; however the model tends to over-develop pressure 
centers. The effects of terrain and non-adiabatic warming 
were in accord with theoretical expectations; however, the 
excessive development dominated the numerical verification 
so that no significant improvement resulted from the in- 


clusion of these refinements. 
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tell RebUCTION 

The availability of high-speed large-capacity elec- 
tronic digital computers has long since made numerical 
weather prediction an operational reality. Notable success 
has been achieved at the 500-mb level. However, surface 
prognoses obtained by numerical integration techniques on 
the average have not equalled the accuracy obtained by 
subjective methods. Even with some relatively complex 
models, the degree of success has not corresponded to the 
degree of sophistication of the models. To improve the 
quality of surface prognosis, two alternatives seem to be 
feasible. One approach is to deviseseven mere elaberare 
mathematical models in order to take into account as many 
pertinent weather parameters as possible. Perhaps new 
parameters which are not just a ree@mbination of the 
familiar physical variables--pressure, temperature and 
hMaMmidity, will have to be introduced and measumed teva high 
degree of accuracy on a synoptic time and space scale. 
Another alternative is to try quasi-empirical prediction 
models with a reasonable dynamical basis. In this paper 
a Simple gquasi-empirical 1000-mb prognostic model is 
derived and testing results are presented. 

Since 500-mb prognoses from the barotropic model or 
LES variations have been smmte successmm@l, a next Logical 
step would be to devise a method to forecast the thickness 
field; and then through the predicted thickness arrive at 


Bie PEegmests oresome other level. “In this experiment, a 


il 





Pochhigue shor pred CcLringmencutnieKness field between 
1000 mb and 500 mb is combined with results of the baro- 
LEOpic model currently in Use by the Fleet Numerical 
Weather Facility. The predicted 500-mb stream field at 
every hour 1S applied in Curae comcecneraroma en exness 
between levels through a thickness=tendency equation. 
After sufficient iterations, a 24-hour 1000-mb prognos ‘¢ 
contour map is obtained. 

From the hydrostatic viewpoint, pressure is simply the 
weight of the air column above a parcateulas level and 
local changes in pressure may be thought of as the inte- 
grated density change in the air column. The density change 
can be described in terms of the temperature field which, 
in turn, can be expressed as a thickness change. The main 
feature of this model, the thickness-tendency equation, 
takes into account the temperature field, horizontal 
motions, and the vertical motions induced by terrain and 
SUbBAeGeeSeELDess. | 

Actually there is no inherent restriction on the choice 
OCfeEneesOU=n> DrediTCELOnN Equation. ~The EChiveknese—vendency 
equation could be combined with any 500-mb model, though 
perhaps requiring some minor modifications in the pro- 
gramming. Obviously, if another 500-mb model is employed 
instead of the particular barotropic model used here, the 
Boa MeeSwOLn LeStiIng May beequtce dittcuenc, and eaeac Empir— 
ical coefficients may not have the same values as obtained 


in this experiment. 





2. Prognostic Model 
The prognostic model is composed of two prediction 
equations. The fiesta es mC strecam—barotropic vor= 
ticity equation, the other is a thickness-tendency equation. 
With the twisting term and vertical advection of vor- 
ticity omitted, the vorticity equation in (x,y,p,t) coor- 


dinates becomes 


<a #V°VG+h)= —Cf+ s)7-V - 7 


Bue GO the fact thateaesmeceemediace levels in the tropo- 
Spitewe between SU0vangsoUGem> sthne horizontal velocity 
divergence is normally relatively small and frequently may 
permeglected, the vomgereteyeequatron may be written anegese 


form 
3o+yvev(ftf)= o , (2) 


Mere geostrophic approximation of (2) with 


er : = L2 


yields the barotropic gquasi~geostrophic model 


vi(®= Te. fv*e+f) 


where J is the Jacobian operator. Equation (3) has had 
considerable success in predicting the height field for 
SIGpiomouewcanie to the velocity divergence inherent in che 
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SeSOS5erOPpnaic Wind, an Gruremmeonmonly referred to as “spuri— 
ous anticyclogenesis" is introduced. To remove this source 
of error, a non-divergent wind is used as follows: 


V.> —Vek Oe 


With these substitutions, equation (2) becomes: 


V(4)+S(¥, Vi¥tt) = 0 , 


A characteristic of the barotropic model as repre- 
sented by (4) is that very long waves, say number 1, 2, and 
3 with respect to a latitude circle, tend to retrogress at 
very high speeds, a phenomenon not observed in nature. To 
Sonmerol this Crror, anwadditicnal cern Nas Seen itnesedilecar 


ime result is the Helmholtz stream—barotropic equation 
oe. Se o 
a ome Gite y J 


mere 2) is the absolute vorticity and At is a constant. 





The first law of thermodynamics and the Poisson's 


equation can be combined in the form 


Q = oo (dt +V-VE +057) . 


Piso LromerolssSon'sS Equation sence ene Caqudmlon Olmoira ce nt 


follows that 


OH = ral . p Cuvier inte (we) 





which may be used in (6) to give 


ag. Cae = Sa ; 
SO ST 


Pamadkily replacing a in thestirmsitStewo eeums om the left 
by means of the hydrostatic equation, and on the right 


miaoughi the equation Cl: state yields 


= ale -V- Vises Tu — a (9) 


Here Fa 1. 98 om 


C oP ap 


memordered to be constant in space and eime- 





isthe stabpmrrty Taceor andes 


The "vertical velocity" may be considered to consist 
of three contributions, the large-scale vertical motion, 
the terrain-induced vertical motion and frictionally-in- 
duced vertical motion. Here these contributions will be 


represented by assuming GQ) to be of the form 


—~™ 


ra 
G) = O),| 1- a UREVY Het I (10) 


where 





The first term on the right side of (10) is the para- 
bolic distribution frequently assumed to be typical of 
large-scale pressure systems [1]. 

The second term of W represents the terrain-induced 
vertical motion. Here the assumption is made that the sur- 
ma@e vertical velocity Wo LS PrEOpOCEETend tOurne enol’ 2zOhear 
wind and the gradient of terrain height at the surface; 
i.e., Wo = kY* VY H,- Since 4), = -g Ip Wy, it follows that 
(J, = -3 f k VY, VY H,. The terrain-generated vertical velo- 
city is further assumed to decrease linearly with pressure} 
thus we arrive at the form in the second term. 

The third term represents vertical motion induced by 
surface friction. Here Wy) is the "vertical velocity" at 
memmoretne friction layer save ror simrace es enesco ems 
assumed to have a simple distribution with respect to 


height which can be approximated linearly as shown in fig. 1. 


() = 30 at “R. 





G) = Cy al p, 


GJ= OC cit ie 


Figure 1:7 Protile ob@vertitcarvelocrtey duermEo 
surface stress. 
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Cressman [2] deduced the following form for GQ), in terms 


of geostrophic surface wind and a drag coefficient: 

y= By (ule )- $(Gulere)|= 42 e an 
where 

—— ay (Cau fos v*)— $§ (ay are] (12) 


Here uand v are the x- and y-components Of the Gecostecpaic 
wind at 1000 mb, f is the coriolis parameter and Ca is the 
drag coefficient. The latter is mainly a function of the 
roughness of the underlying surface and stability in the 
micrometeorological layer, and has been determined for 
average conditions by Cressman for large-scale forecasting 
evem the Northern HemisSpnere. 

Upon integration between ie and ile equation (9) 


becomes 
> | 3 _- 

Sh Vo Vh- FO hes - ERIC AVE THe 
= + 0 7 IE. +06tcA q Q 


where Q represents Ene aMoumt gOruneat cOdCd See mtincm layer 


(13) 


per unit time. 
If in equation (13) we 


1. assume Of RV-Vh and combine it with 


the first term and affix a coefficient K); 


\ 
‘. 





2. combine coefficients in the second term and 
designate the combination as Kyi 

3. assume Q is igre Posies eieul to Wey T Ol 
0 = KV eV Ts: where T. is the sea-surface 
CEMpeEauuTe wane 

4. assume a reasonable pressure level as top 
of friction layer in the frictional term 
and designate the combined coefficient as 
ie 

we arrive at: 


— -KVs Vk -KGVeV HE + V-V Ke 04 


tis its the thickness—téendéency equation, witherine Ccmeimareal 
Goelficients to be determined through testing, Theventc.— 
ness equation together with the 500-mb stream-barotropic 
equation constitute the prognostic model. 

The local rate of thickness change is thus seen to be 
due to the combined effect of thickness advection between 
layers, terrain effects, non-adiabatic heating and a 
frictional effect. Since cold advection is usually 
associated with adiabatic warming due to subsidence aloft 
and by heating from below in the surface layers, while 


warm advection is usually counteracted by cooling due to 


large-scale upward motion and in the surface layers by 


Gooling from below, the coefficient K) might be expected 
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to have different values for cold and warm advection. This 
was borne out by the results of this numerical experiment. 

The terrain term represents vertical motion induced 
by uneven terrain, and obviously it vanishes over oceans 
and level land. 

The non-adiabatic term in the form KWVeVT. represents 
vertical diffusion of heat between air and surface water. 
For an air maSS moving over warmer water, the surface 
temperature quickly adjusts to that of the surface water 
and rapid vertical diffusion of heat takes place. The rate 
of diffusion will normally increase with increasing hori- 
zontal wind which enhances the turbulence, and upward dif- 
fusion of heat will certainly continue as long as the water 
temperature increases downwind. For warm air moving over 
colder water the stability created by surface cooling in 
the lowest layer will tend to inhibit diffusive heat ex- 
change. This suggests the empirical coefficient K. may 
also be assigned two values depending on the sign of 
VV T,- 

The effects of surface stress about closed pressure 
systems may be summarized as follows: 

Low: upward motion at top of friction layer, Wy 
negative; conEribution to thickness tendency, 
negative. 

High: downward motion at top of friction layer Q)y 
positive; contribution to thickness tendency, 


Posiolve. 





Hesselberg and Fricdmann’ found=that thewfirietion efrece 
is at least two orders of magnitude smaller than some of the 
Gther terms in the vorticity Cquaeren a .eeosSunpti One. 
made that it will also be quite small as compared to other 
terms in the thickness-tendency equation. Hence in the 
present experiment it was neglected. Thus equation (14) 


= ~ RV VR-RVE VE TKVe'V Ty (15) 


where fs represents the wind at the terrain height. Next 


the winds V. and Vy) in equation (15) are approximated by 


5 


geostrophic values resulting in 


ab = - $ KT, AK TW] KM-VH 08 


The wind at terrain height is computed by linear inter- 
polation. On the assumption that wind hodograph is approx- 
imately a straight line between 1000 and 500 mb, V, may be 


expressed in the form 


V=V- (A444) A 
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Thus the wind at terrain height is equal to the wind at 
500-mb level minus some fraction of the thermal wind, and 
the fraction is the ratio of the surface/500-mb thickness 
to the 500/1000-mb thickness. With the geostrophic approx- 
imation for the wind and the letter A to designate (2s-He | 


Vv 
Hee LtOllows that 








Vi=+(- 48i+ 2255) -A-—-(- 2a 4 2k 
cay Jy aN ak 5 Pat as / (18) 

Hence the terrain term in equation (16) may be computed 

using two Jacobians. However, because of limited computer 


storage capacity, instead of the two Jacobian operations, the 


following transformation is made 





Veh HOSS S)-AC SHE DPVH, a9 


In finite differance form, the last expression becomes 


Ve VHF {Lo (SurAhad (Zoo A hd(Hee-Hew) 


+ Ze Ahp) ( Zsw" had (Hy Hes) (2 


0) 


il 





where the subscripts N, S, E and W are notations for the 
location of grid points. The quantity A has three var- 
iables, but may be simplified by taking standard atmos- 
pheric values for 2, and h without introducing too large 

an error. With the letters Zine Zig: Zee and Ziw to denote 


the quantities within the brackets in the final expression 


for VV, above, equation (15) has the form 


pA 


at apa (KT (2p WOKS (2.7% 


® 


a2 





3. Procedure 

Three days’ data invapmeel Sewers Lested. Progne— 
ses were made for April 2nd, 3rd and 4th and were labeled 
as Cases 1, 2 and 3, respectively A CDC-1604 digital 
Somputer was employed fer ene Conpmeartonssover a 1977 
point octagonal grid for the Northern Hemisphere poleward 
of 1ON. At the end of each iteration, a predicted thick- 
ness is obtained at each grid point. This is Subtrmacced 
from the prognostic 500-mb stream function to find the 
1000-mb height. The final output is a predicted 1000-mb 


height field in a contour map form. The forward-difference 





method was used for calculation of Dh for the first hour 


a 


and a centered-difference method was used for succeeding 


hours. The forward jand centered dit ference methods may ase 


expressed respectively as 


i. (2 At 


fee sat Wisk] HE 


and 


I/F 
bea ye = heat t (SH) 240 


A schematic block diagram of the thickness prediction 
PEOgGgeEam 1S Given in figaaze 


Using the computer, a "pillow" and a RMSE for each 
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PRESET 


GENERATE gm 
Afa2 
READ IN Ze | 
COMPUTE : 
Kj J(Z. ,h) 


Ko Vv, - VY AL 
— J(Zo, We 


COMPUTE 2) py 
Bye 


adding 3 terms 
and emul eile ing 
gm 3 and 3600 sec. 
afd 





GENERATE h 


from 2h_ 
ai 





GENERATE 
HOURLY Z,. 
( Z 

O 





PRINT MAP 


Figure 2. Schematic block diagram showing procedure. 
14 





prognostic map were computed by compdrison with the actual 
map for verification purposes. The formulae used in com- 


puting the pillow and RMSE are 


xX 
= (A-B),, 


Fe 


pillow Lite: 


S$ [(4-8)- pi (Loup 


X 








RMSE = 


where A and B stand for predicted and verifying values, 
respectively, at the same grid poeintweand x stands for) the 
total number of points differenced. Qualitatively, a pillow 
is the algebraic mean difference over the whole grid between 
the two fields A and B; while the RMSE is a numerical meas- 
ure of the prognosis, and it reflects errors in position 

and in intensity of pressure systems. 

Due te fictitious retlectiememenstrne boundary pemnes 
and inaccuracies in the geostrophic approximation at low 
latitudes, it was decided to limit the verification to the 
areacimorsenesof tiv 20th Jlatreudeverrecle- 

Since the major contribution to thickness change is 
due to the advection term, it was investigated first. 


Values of K, determined by previous investigators [3] were 


IS 





used initially, then the value was varied to arrive at the 
best verification. Furthermore, Since warm and cold ad- 
vection did not necessarily give the same contribution 
EemrEne predicted change in thickness, combinations of eye 
different values were tried during the investigation of 
Ears term. 

With the coefficients for the advection term tenta- 
tively determined, the investigation was carried on to 
find suitable coefficients for the terrain term. Only one 
coefficient of best fit was sought for the terrain term. 
MPHring the investigation of this term, a smoothed terrain— 
netgm:e Lield 1S required IM order erouoerainschie Signi. 
cant scale of vertical motion. Since the grid distance is 
about 380 km, any smaller-scale variation in the terrain 
field would not be included. For this reason, the smoothed 
terrain-height field data of J. Smagorinsky was used. 

Finally, an investigation was conducted to find suit- 
able coefficients for the non-adiabatic term. Two coef- 
ficients were sought for this term, one for heating and 
one for cooling. Values of K., determined by previous inves- 
tigators [4] were also used initially in this case. The 
monthly-mean Sea-temperature field of April was used during 
the investigation of the non-adiabatic term. This field 
was obtained by graphical addition of the monthly-mean 
air-temperature charts of April and the 3-month-mean air- 


sea temperature difference charts [5]. The resulting 
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temperature fiela 1s Shown inwesg. ve. The temperature 
field on the land area was included to facilitate the 
computations; however, these values did not enter into the 


final results. 


1 





4. Results and DiscuSSion: 

Figures 3 to 8 indicate the initial maps and the 
prognostic maps using different combinations of the terms 
forming the sthickness-tendency Cauatelcneaowtsolines 
are drawn at intervals of 200 feet, and labeled in hun- 
dreds of feet. Case 3 has been chosen for illustration 
purposes because it gave the best results among the three 
cases tested, both in the prediction of the positions of 
Emempressure centers and an numenmiecaleverriieo4-toemcr 

By comparing figs. 4 amd See it can be Scene rhaceeic 
positions of the pressure systems are in close agreement 
with four exceptions, namely: the high over the east coast 
of Asia; the high off the west coast of the United States; 
mae high over the Caspian ™Séa;"andy the higheilying toed 
west of Gibraltar. Possible reasons for errors in their 
locations will be discussed at the end of this section. 

It was also of interest to notice that the positions of 
systems in the Atlantic Ocean were predicted more accurately 
than those in the Pacific Ocean. 

Before discussing the coefficients found thus far in 
this investigation, it iS necessary to mention two phenom- 
ena related to development. These phenomena are fictitious 
anticyclogenesis and cyclogenesis which predominated in 
all three cases tested. While no dynamical explanation 
could be determined from this investigation, it was apparent 
that the effect was due mainly to the advection term. 


Comparison of figs. 4 and 6 shows that most of the highs 
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were too high by an average amount of 200 feet, while most 
of the lows were too low by an average amount of 100 feet 
when the advection term was used alone to make the forecast. 
since the major contribution to thickness change is due to 
the advection term, these phenomena Pavarrtawinvee concealed 
the effects of a refinement like the non-adiabatic heat 
exchange in the numerical verification. 

Aside from the effect of over-prediction by the advec- 
ELON term, it was found that the combination of cocrr. 
cients 0.4 for warm advection and 0.5 for cold advection 
gave the best verification when this term is used alone to 
make the forecast. Other sets of coefficients, such as 
O076/70.9, 0.7/0.8, 0.5/0.5, 0.3/0.4, tormention only ate. 
gave larger RMSE values. Therefore, tegwas deci dedsrenca. 
tively to use 0.4/0.5 as the coefficients of best fit for 
the advection term. The results of numerical verification 
BOG@EpEOGNOSeS uSing the advection Temnmyonly wien cecrrta. 
cients of 0.4 for warm advection and 0.5 for cold advection 


are shown in Table l. 
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Case Pillow RMSE 


Number (ft) ) (Ee 


1 
(0000Z April 2, | ° : + 189 
1955 
: 
De ie 
(0000Z April 3, iy | + 
1955 ) 
: 
: : 
3 
ae ) + 168 


(0000Z April 4, 
1955 ) 





Table 1. Results of numerical verif™@®€ation for prognoses 
using the thicknegs advection term only with 
coefficients of OU.) for warm @e@vection and 0.5 
for cold advectign, 
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For the investigation of the terrain term, verifica- 
tion by pillow and RMSE in Table 2 showed no improvement 
in RMSE when the terrain term was included in this model. 

Im this connection, it should be noted €hat the empirically-— 
determined best-fit coefficient of the terrain term is 0.01, 
while the maximum terrain-height gradient is about 12 times 
larger than the thickness gradient. Hence the terrain term 
contribution can be at most about 25% of the advection term. 
Also, less than half of the total area covered by the grid 
point is affected by the terrain term, while the pillow 

and RMSE were computed over the entire grid. Therefore, 

the actual contribution of the terrain term should be 

larger than that shown by RMSE. 

In order to investigate the effects of the terrain- 
induced vertical motion in more detail, a differential 
analysis was made and presented in fig. 9. Here the dif- 
ference between the forecast with thickness advection only 
(fig. 6) and that with both the terrain term and thickness 
advection (fig. 7) was computed. The isolines are drawn 
at intervals of 50 feet, and the figures labeled in 
tens of feet. Of the 20 closed isolines observed in fig. 9, 
14 indicate a decreased error; 5 of them an increased er- 
ror; and one, neutral. Thus, there are some areas where 
the terrain term does not give any improvement and, in 
fact, made the results worse. For example, over Greenland, 
the forecast made with the thickness advection alone over- 


intensifies the low, and with the inclusion of the terrain 


2 





SD vere ep ap es oe - 


(OO00Z April 2, 


Case 


Number 


955 


(OO00Z April 3, 


(OCOOZ April 4, 
1955 


Table 2. 


2 


1S) 
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term, further deepening is observed on the leeward side of 
the mountain, giving an even poorer result. 

Fourteen closed isolines of the 20 on the differential 
analysis showed a better forecast, with a two-to eight-mb 
improvement in the surface pressure field, when the 
terrain-induced vertical motion over land areas was in- 
cluded. These results indicate some justification for the 
inclusion of the terrain term in the model. 

Since the terrain-induced vertical motion is a) fune= 
tion of the wind component parallel to the terrain-height 
gradient as well as the terrain-height gradient itself, 
large values of terrain contribution may be expected in 
areas where above-mentioned values are large. In this 
respect, the test data available for this investigation 

were not favorable for a good evaluation of the terrain 
term. The synoptic maps show that the areas where the 
terrain gradient is large are mainly associated with 
siciilewhind Components paballel jo the terra n-nerge 
gradient. 

With only three days' data tested, no final conclu- 
sion can be reached regarding the terrain-induced verti- 
cal motion in this model at this time. Effects of the ter- 
rain-induced vertical motion were produced in areas where 
they were expected from dynamic considerations; however, the 


values are too small in general. This seems to suggest 


PA) 





that further experiments may indicate a larger coeffi- 
cient for the terrain term with perhaps some other addi- 
wive empirical factors. 

The results of numerical verification for prognoses 
using advection and non-adiabatic terms are shown in Table 
3. The eight blanks in Table 3 arose in cases where the 
forecast map was very much distorted, in contrast to the 
regular appearance of the remaining forecast maps. This 
erratic behavior was apparently connected with a malfunc- 
Elomee: Che magnetic tape unres- 

mae coerricientsS Of best fie for Cheowmen-cclasaere 
term, as shown in Table 4, were 10 for heating and 6 for 
cooling. They were chosen only because they provided the 
smallest pillow. The RMSE obtained by using this set of 
coefficients was +168 feet for case 3, which equaled 
exactly that obtained by using advection term alone for 
prognosis. As a matter of fact, Table 3 shows that the 
RMSE values increased when the coefficients were increased 
algebraically. Since it was mentioned earlier that the 
RMSE is a numerical measure of the accuracy of the forecast, 
fPeeappears that the addition Ob Ene mon-adwabatic tern ac 
best made no improvement in the forecast. Two arguments 
meeceto be made clear at this point. )FPirstiy, fictitious 
amntteyclogenesis and cyclogenesis by the advection term 
alone already have over-intensified most of the highs and 
the lows. Since the effect of the non-adiabatic term on 


the average is one order of magnitude less than that of 
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Table 3. Results of numerical verification for prognoses 


using advection and non-adiabatic terms with 

coefficients of 0.4 for warm advection and 0.5 
for cold advection. The coefficients used for 
the non-adiabatic term are shown in the table. 
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Table 4. Results of numerical verification for prog- 
noses using all three terms of the thici:ness- 
tendency equation, with coefficients as shown 
in the table. 
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the advection term, the advantage of including the non- 
adiabatic term would hardly appear. Secondly, for the non- 
adiabatic heat exchange, the zones of thermal concentration, 
where transformation of polar-continental air to polar- 
maritime air takes place, are really of primary concern. 
During the month of April, these zones of thermal concentra- 
tion were found along the northeast coasts of Asia and North 
America. For the three cases investigated, cyclonic activ- 
ity prevailed in these areas. With the 1000-mb lows al- 
ready over-intensified by the advection term and with the 
flow oriented such that heating took place in these areas, 
the lows were intensified further. For example, the low 

off the east coast of Japan in case 3, where the lowest 
observed height was +30 feet, the prognosis with the ad- 
vection term alone gave a lowest height of - 170 feet while 
Enewnon—-adiabatic term decreased “ene lowest valueminrieaer 
EOw= 9250 feet. Thus, in Spite of the face enat ene bens 
i10r was in accord with theoretical expectations, the RMSE 
Gree netesimprowe when the non-adiabatic term was included 
ines case. 

In order to investigate the effect of the non-adiabatic 
term further, a differential analysis was made by sub- 
tracting a prognostic map with both advection and non- 
adiabatic terms included (fig. 8) from the one with 
Only the advection term.) This 1selaceq ene vert ece erence 
non-adiabatic term (fig. 6) and it was found that by using 


a coefficient of 10 for heating, the largest 24—-hr 


| ied 
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contributions by this term in areas of thermal concentra- 
tion were -282 feet off the east coast of Nova Scotia, at 
41N, 55W; and -230 feet off the west coast of Japan, at 40N, 
134wW. 

Two possible sources of error warrant discussion. In 
the first place, the 1000-mb map was a relatively flat map, 
with the range of thickness values and 500-mb height much 
greater than the range of the 1000-mb height. Therefore, 
errors in the 500-mb height, when transmitted downward 
through the thickness, can cause serious deviations at the 
lower surface. Furthermore, it was obvious that one could 
not expect the locations of the pressure systems on the 
lower surface to be better than what the 500-mb prognosis 
can provide for this model. Secondly, the fact that the 
POsiticnms, Of the pressure SyStems were ai bert rer ee need 
Atlantic Ocean than those in the Pacieice Ocean, t20c. mime 
three cases, was most probably due to the fact that weather 
reports in the Atlantic Ocean were more numerous and accu- 
rate than those in the Pacific Ocean. In other words, the 
initial data from which the prognoSis was started could 


have been in error to a certain degree in the Pacific Ocean. 





Se cone lust ons: 

starting with a dynamical basis and utilizing the 
rapid computing capabilities of the CDC-1604 computer, 
empirical coefficients were found for each of the three 
terms of the prognostic equation designed to predict the 
1000/500-mb thickness. By subtracting the predicted 
thickness from the 500-mb 24-hr prognosticated stream 
function, a 24-hr prognosis of the 1000-mb map was obtained. 
the fact that this model did provide avmap wiLenemoseyors en. 
Pressure systems appearing at ERG Eaghesplaces 1S Cneeumeag— 
ing. Due to the relatively small Sample tested, the empir- 
ical coefficients found thus £ar were not necessarily the 
best ones; but they will nevertheless provide a starting 
point for further investigations of a Similar nature. 

Improvements of the model might include the incorpor- 
ation of certain empirical rules in order to remedy the 
phenomena of over-predicting cyclogensis and anticyclo- 
genesis and the use of the weekly-mean sea-temperature 
field instead of the monthly mean. Furthermore, any im- 
provement in the prediction of the 500-mb maps will give 
increased accuracy at 1000 mb. In conclusion, it is the 
opinion of the present investigators that this model can 
provide a first approximation in the preparation of the 


surface prognostic map. 
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figure 3. The 1000-mb map for 0000Z April 3, 1955. 
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Figure 6. The 24-hr prognostic map with the thickness 
advection term only. 
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Figure 7. The 24-hr prognostic map with the 
thickness advection and the terrain 


terms. 
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Figure 8. The 24-hr prognostic map with the thick- 
ness advection and the non-adiabatic 
terms. 
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Figure 10. The monthly-mean sea-temperature field 
GE April. 
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